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Abstract. Two Delone sets are called homo-metric when they share the same autocorrela- 
tion or Patterson measure. A model set A within a given cut and project scheme is a Delone 
set that is defined through a window W in internal space. The autocorrelation measure of 
A is a pure point measure whose coefficients can be calculated via the so-called covariogram 
of W. Two windows with the same covariogram thus result in homometric model sets. On 
the other hand, the inverse problem of determining A from its diffraction image ultimately 
amounts to reconstructing W from its covariogram. This is also known as Matheron's co- 
variogram problem. It is well studied in convex geometry, where certain uniqueness results 
have been obtained in recent years. However, for non-convex windows, uniqueness fails in a 
relevant way, so that interesting applications to the homometry problem emerge. We discuss 
this in a simple setting and show a planar example of distinct homometric model sets. 



1. Introduction 

Quasicrystals form an interesting class of solids whose structure determination requires 
a substantially extended setting in comparison to classical crystallography, see j^l EH an< ^ 
references therein. In particular, one first has to extract from the diffraction data the higher- 
dimensional embedding space and an appropriate lattice in it, followed by the reconstruction 
of the window (or acceptance domain) for the description of the atomic positions. 

In what follows, we begin to analyze the structure of this inverse problem from a math- 
ematical perspective. To this end, we assume a perfect diffraction experiment and the ap- 
propriateness of describing the underlying structure as a model set. The first step, finding 
the Fourier module and then the resulting cut and project scheme, seems relatively straight- 
forward. We thus assume this step already done, see also jS] and further remarks below for 
more. 

Consequently, we need to determine the window to complete the picture. Here, we discuss 
the question how and to what extent it is specified by the diffraction image. Since all princi- 
pal aspects show up already for a mono-atomic system, we make this simplifying assumption, 
too. We are well aware of the fact that real quasicrystals require the reconstruction of several 
windows, but this adds a technical complication, not a principal hurdle, to the picture de- 
scribed below. Complementary reconstruction methods on the basis of discrete tomography 
look also promising in the context of quasicrystals 4 , but are not considered here. 

A Delone set A is described by the corresponding Dirac comb 8jy := J2 x eA^, where 5 X 
is the normalized point (or Dirac) measure at x, see for details. A perfect diffraction 
image of A, as described by the positive measure j A , uniquely determines its inverse Fourier 
transform, which is the autocorrelation (or Patterson) measure j A . Our starting point is 
thus the (hypothetically complete) knowledge of 7^, where we shall work with infinite point 
sets A to avoid the extra layer of complication from finite apertures. The remaining task is 
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then to determine a window W from this information, which leads to Matheron's covariogram 
problem 0. 

This is an example of a class of inverse problems where one aims at the reconstruction 
of a finite or compact set in Euclidean space from its (possibly weighted) difference set. As 
originally observed by Patterson |18j for finite point sets in M, there need not be a unique 
solution to the reconstruction problem. To capture this ambiguity, two point sets are called 
homometric when they share the same (weighted) difference set, see |22) and references therein 
for further examples. 

Similar questions emerge for non-empty compact subsets K C M. d , where the covariogram 
g K {x) = vol(-KT n (x + if)) encapsulates the difference information, and one tries to determine 
K from the knowledge of g K , see jZlE]- Furthermore, one is also interested in similar concepts 
for infinite point sets in general, and Delone sets in particular. Here, two Delone sets are 
called homometric when they share the same autocorrelation measure, to be defined in detail 
below. There is an interesting connection between the homometry of Delone sets and the 
covariogram problem for the class of model sets (also known as cut and project sets) |16j . 
which is explored below. 

The paper is organised as follows. The covariogram and its elementary properties are 
introduced in the next section, while Section|2]recalls the basic facts about the autocorrelation 
of a regular model set and links it with the covariogram. Within a given cut and project 
scheme, the homometry problem thus becomes equivalent to Matheron's covariogram problem 
(Zj. Section 0] shows a planar example of two different sets with the same covariogram, which 
give rise to two distinct homometric model sets. 

2. Properties of the covariogram 

Let K C M. d be a non-empty compact set which is the closure of its interior, K = K° . The 
function 

(1) g K {x) := vo\(Kn(x + K)) 

is called the covariogram of the set K. This defines a real- valued function on R d which is 
inversion symmetric, g K {—x) = g K (x), and satisfies 

sup g K (x) = g K (0) = vol(if). 

The function g K {x) is continuous on M d , see ^3 Thm. 3.1]. Moreover, g K is a positive definite 
function with compact support. The latter is the so-called difference body 

(2) suppfe) = K-K := {x - y \ x,y G K} , 

which is always inversion symmetric. In fact, the entire covariogram of K equals that of —K 
as well as that of any translate K + t. This means that g K can determine K at best up to 
translations and inversion. As we shall see below, even this is generally not the case. 

If 1 K denotes the characteristic function of K, the function g K (x) is given by the convolu- 
tion 



(3) 



9k( x ) = (J-x* 1 -^)^) 
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Figure 1. Example of two homometric polyominoes, each with an equally 
shifted copy overlayed. This illustrates that and how the two intersection 
areas coincide. For their use as windows in the cut and project scheme, we 
choose the marked dot as the origin in internal space. 



The Fourier transform thus satisfies 

(4) 9k (k) 



which is a positive, analytic function that vanishes at oo. This relation is the reason why, if 
K is itself inversion symmetric in the sense that —K = t + K for a suitable translation t, 1 K 
can be reconstructed from the knowledge of g K , up to translation and inversion |1U| . 

If K is a convex polytope in dimension d < 3, it is determined by g K , in the sense mentioned 
above, see [SJ El IE] and references therein. Other examples include convex sets in the plane 
with C 2 -smooth boundaries, and large other classes (e.g., non-C 1 bodies and bodies that are 
not strictly convex), though the general claim is still open, see [7j. However, starting with 
dimension 4, uniqueness fails even within the class of convex polytopes, compare [Tj^j. In 
general, the reconstruction of K from the knowledge of g K is a difficult problem, which is 
beyond our scope here. However, an interesting example of two polyominoes with the same 
covariogram was constructed in ^2] an d is shown in Figure ^ 

3. HOMOMETRY OF MODEL SETS 

In this section, we use the setting of Euclidean spaces for simplicity. This covers the most 
common situation in practice; we refer the reader to |16[ I20| for the more general setting. 

The basic cut and project scheme ^U] requires a number of spaces, sets and mappings as 
follows. 



(5) 
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L 



1-1 



U lattice 

c 



(J dense 
L* 



L 



4 



MICHAEL BAAKE AND UWE GRIMM 



Here, M. n and M. m are called the direct (or physical) space and the internal space, respectively. 
The mappings ir and 7r int denote the canonical projections. The set £ is a point lattice in 
R n xR m = M n+m such that n is 1 - 1 between C and L = vr(£) and that L* = vr int (£) is dense 
in R m . Finally, the *-map is defined on L via t* = ^int (C 7r lz:) 1 (*)) f° r an t E L. It has a 
unique extension to the rational span QL of L, but not to all of M. n . 

To define a (regular) model set on the basis of the cut and project scheme ©, we need to 
specify a window W C M. m . We assume that W is a compact set with non-empty interior and 
the property that W = W°. Regularity refers to the extra requirement that the boundary 
dW has measure in internal space. In this setting, we obtain a (regular) model set A as 

(6) A = t + X(W), 

with any t £ R n and X(W) := {x £ L \ x* G W}. 

Let us now consider the autocorrelation 7^ of the Dirac comb 5 a = YlxeA^^ which is 
defined as 

(7) 7 A ■= lim 7rB^ S ^B r *S- AnBr , 

r->oo vol(i3 r ) 

where B r denotes the open ball of radius r around 0. The limit exists by general properties 
of model sets, compare |2U| . Moreover, it is independent of the translation t in ©, wherefore 
we set t = from now on without loss of generality. It is well known that the autocorrelation 
has the explicit form 

(8) 7a ■= Vi^t* 

xeA-A 

where A— A is a uniformly discrete (and hence countable) subset of M. n and the autocorrelation 
coefficient r](x) is given by ^1 Prop. 3] 

(9) v(x) = dens(yl) v ° l{w 2m)~ x)) 

= dens(£) vol(W n (W-x*)). 

This expresses the autocorrelation of A in terms of the covariogram of the window W, see 
[21 E] for related problems that can be expressed this way. 

Let us recall that two Delone sets are called homometric when they share the same auto- 
correlation, see for an example in the setting of quasicrystals. Note that this definition, 
due to the infinite volume limit in Eq. (|7|). disregards contributions from point sets of den- 
sity 0. For simplicity, we restrict ourselves to Delone sets with unique autocorrelations here 
(meaning independence from the selected averaging van Hove sequence, compare [SOI)- In 
particular, this is the situation for all model sets. 

Theorem 1. Two model sets obtained from the same cut and project scheme are homometric 
if and only if the defining windows share the same covariogram. 

Proof. If the two model sets, A and A', possess the same autocorrelation, the coefficient rj(x) 
is valid for both of them. Since the cut and project scheme is fixed, the factor dens(£) in 
Eq. © is the same for A and A' . Consequently, the covariograms of the windows coincide on 
all points x G A — A. By the general setting of a cut and project scheme, A* — A* is a dense 
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subset of W — W, so that the covariogram, which is a continuous function ^U], is completely 
determined on W — W, while it vanishes on the complement. As the same argument applies 
to W' — W, the covariograms of W and W coincide. 

The converse claim is obvious from Eq. ©. □ 

It is a well established result, see |20| and references therein, that the diffraction measure of 
a regular model set A is a pure point measure, meaning that it is a countable sum of weighted 
Bragg (or Dirac) peaks. It is given by Ja, the Fourier transform of the autocorrelation measure 
7*. Since both 7^ and 7^ are tempered measures, and since Fourier transform is invertible 
on this class of measures, the following consequence is immediate. 

Corollary 1. Homometric model sets from the same cut and project scheme have the same 
diffraction measure. Conversely, kinematic diffraction cannot discriminate between homomet- 
ric model sets. □ 

So far, the formulation was based upon a fixed cut and project scheme. In general, however, 
there are different possibilities to generate a given model set. This also implies that two regular 
model sets obtained from different cut and project schemes can still be homometric. If so, 
they have the same autocorrelation 7 by definition, and thus the same diffraction measure 7. 
Since regular model sets are pure point diffractive ^UEO], 7 comprises a countable sum of 
Dirac measures with positive intensity. The Z-span of their positions is the so-called Fourier 
module. Via minimal embedding, it defines an essentially unique cut and project scheme, 
whose dual in the sense of ^H] is a natural setting to describe both model sets in the same 
cut and project scheme, see [3] for details. In this sense, the restriction in Theorem ^ is not 
essential. 



4. A SIMPLE PLANAR EXAMPLE 



A simple example of two distinct homometric polygons was given in ^21- They are shown 
in Figure together with some typical intersection pattern. The polygons Pi and P2 are 
polyominoes, which both consist of 15 elementary squares. Their joint covariogram, supported 
on the inversion symmetric difference body, is sketched in Figure 121 as a contour plot. 

To turn this pair of polyominoes into a pair of homometric model sets, we use them as 
windows for L = Z[£] C M 2 with £ = exp(2-7ri/8), with the internal space M 2 and a *-map 
defined by a suitable algebraic conjugation ^2]- Here, we use the one induced by £ 1— ► 
With 

this choice amounts to the replacement y/2 1— > —y/2 on the level of the coordinates. 
The model set formulation is thus based on a lattice £cR 4 . Its basis matrix reads 



(10) 



B, 



f 1 1/V2 





-1/V2 \ 


1/V2 


1 


1/V2 


1 -1/V2 





1/V2 


V -1/V2 


1 


-1/V2 J 



where the columns contain the coordinates of the basis vectors. The first two rows show the 
coordinates in direct space, the remaining rows those in internal space. Since |det(5 £ )| = 4, 
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Figure 2 . Contour plot of the covariogram of both Pi and P 2 ■ It is inversion 
symmetric with respect to the origin, the latter marked by a dot. The outer 
polygonal line is the boundary of the difference body, which is the support of 
the covariogram, see Eq. (J2J). 

the lattice C has density dens(£) = 1/4. In this setting, we define the two model sets 
(11) Ax = A (Pi) and A 2 = X(P 2 ). 

Both windows are placed in internal space such that their lower left corner has coordinates 
(—1/2, —1/2). This way, the corresponding model sets are generic, meaning that L*PidPj = 
for j G {1,2}. They have density 

dens(yli) = dens^) = ^ 

and are homometric by construction, due to Theorem Note that Ai and A 2 are not in the 
same LI class, compare and that they differ in points of positive density, see below. 

A patch of A\ is shown in Figure 01 while Figure 0] displays the difference between the two 
model sets. This can be understood in terms of the windows, compare Figure[5J In particular, 
the difference sets are model sets themselves, with 

dens(A 2 \Ai) = dens(Ai\A 2 ) = ^dens(vli) = \. 

Note, however, that the sets A 2 \ A\ and yli \ A 2 are not homometric. 

The lattice dual to C is denoted by C* (observe the different star) and is spanned by the 
basis matrix 

B c . = {B^Y = \B C . 
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Figure 3. Circular patch of the model set A\. 



• o • • 



o 



Figure 4. Comparison of the two patterns A\ and A2. Big dots mark points 
that are common to both, while open circles (small dots) are points of A2 \ A\ 
(ofA 1 \A 2 ). 
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Figure 5. Comparison of the two windows P\ (grey area) and (bounded 
by black line). The two white squares on the left (the two grey ones on the 
right) form the window for the points in A2 \ A\ (in A\ \ A2). 

Its projection, tt(C*) = ^L, is then the Fourier module that carries the Bragg peaks. The 
diffraction measure 7 is the same for both A\ and A2, and reads [TT 1141 |2*U] 



E w 5k ■ 



The intensity function is given by 



I(x) = (dens(A,-)) : 



a 2wixy 



dy 



VOl(Pj) J Pj 

where j S {1, 2}, with I(x) being the same function for both choices. Recall 



Jlmrs 



ds 



,7vir(a+b) sin(7rr(6 - a)) 



irr 



the linearity of Fourier transform, and its factorisation over rectangular domains. With 
k* = (k, A), one obtains 

2 



I( K ,X) 



f(n,\) /sin(7TK) sin(-7rA) 



16 



tt 2 k\ 



where 



/(as, A) = (3 + 2 cos(2vrA) + 4 cos(vrA) cos(tt(2k + 3A))) 
(5 + 6cos(27tk) + 2cos(47tk) + 
4 (2 cos(7tk) + cos(37tk)) cos(-7t(3k + 6A))) . 
The diffraction image is shown in Figure El 
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Figure 6. Circular portion of the diffraction pattern, valid for both A% and 
A2- A Bragg peak at k is represented by a disk centred at k with an area 
proportional to the intensity. Note that the image is inversion symmetric, but 
has no reflection symmetry in a line. 



5. Discussion and outlook 

The characterisation of homometric model sets is possible via the window covariogram. 
However, the corresponding homometry classes do not consist of model sets only, since one 
can always add and subtract point sets of density 0. More interestingly, there is still the 
possibility for other sets in these classes that differ from standard model sets in an essential 
way. A better understanding of the homometry classes defined by model sets and related 
structures with pure point diffraction thus requires further investigation. 

When considering homometry for systems with mixed spectra, in particular those including 
diffuse scattering, the situation gets significantly more involved. As an example from |13|. we 
mention the binary Rudin-Shapiro chain in one dimension, which is deterministic, versus the 
binary Bernoulli chain, which is stochastic. They define a homometric pair of point sets that 
even differ in entropy, being versus log(2) in this example. Consequently, new concepts for 
further progress are required. 
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